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ABSTRACT 


The semi-Lagrangian, semi-implicit method 1s used to model the one dimensional shal- 
low water system of equations with surface topography. The forecasts are compared to 
finite difference and semi-Lagrangian, explicit forecasts. In the first experiment, a non- 
rotating system is considered. The semi-Lagrangian, semi-implicit model agrees very 
Well with hydraulic jump theory, while the semi-Lagrangian, explicit model exhibits ex- 
cessive smoothing and the fimite difference model breaks down when the nonlinear 
interactions become too large. In the second experiment, the system is allowed to rotate 
to examine the effect of rotation on the formation of topographicallv induced hvdraulic 
jumps. Although further study is necessary, it 1s clear that rotation retards the devel- 
opment of the low pressure to the lee of the obstacle. A larger domain and higher spaual 


resolution are needed for more detailed simulation of hydraulic jumps. 
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I. INTRODUCTION 


There are many numerical methods available for atmospheric modelling, but two 
features of the semi-Lagrangian, semi-implicit (SLSI) method seem to indicate that it is 
particularly well suited for numerical weather prediction. 

The first of these traits is efficiency. Because the semi-Lagrangian method approxi- 
mates the total time derivative and the terms responsible for the high frequency gravitv 
waves are treated implicitly, the time step is not limited b_ the Courant-Friedrichs-Lewy 
stability criterion. Thus, a relatively large time step can be used with the SLSI scheme 
to produce a timely forecast. 

Secondly, a recent investigation by Kuo and Wilhams (1989) has shown that while 
more traditional forecast models behave poorly when large gradients develop in a rela- 
uvely narrow zone, the semi-Lagrangian technique can still produce relatively accurate 
solutions when a scale-collapse occurs. This is an important attribute for manv 
meteorological applications, such as predicting the formation of squall lines and fronts, 
or investigating mountain effects. 

In this study, the SLSI method is used to model the shallow water system of 
equations with topography in order to evaluate the handling of mountain effects by the 
SLSI method. One set of experiments is conducted for the non-rotating system. Several 
sets of initial conditions are used, some of which are predicted to produce hydraulic 
jumps. The SLSI forecasts are compared with the forecasts of two other models: a 
leapfrog finite difference scheme (FDEX), and a semi-Lagrangian. explicit (SLEX) 
model. Next, the experiments are repeated for a rotating system to examine the effects 
of rotation on the formation of orographically produced jumps. It is Known that rota- 
tion can prevent the development of hydraulic jumps in a flat bottom system (Williams 
and Hori, 1970). 

In the next chapter, the development of the semi-Lagrangian technique 1s reviewed, 
along with theory on the formation of hydraulic jumps. In Chapter III, the governing 
equations and formulation of the FDEX, SLEX, and SLSI models are discussed. The 
initial conditions and results of the experiments are presented in Chapter IV, and con- 


clusions and recommendation for further study are in Chapter V. 


HW. BACKGROUND 


A. THE SEMI-LAGRANGIAN SCHEME 
In a pure Lagrangian scheme, such as Fjortoft (1952) proposed, one set of fluid el- 
ements is tracked for the entire integration period. Thus, a parameter, Q, that is con- 

served following the fluid such that 
dQ 60 


— + V 
Pe acs 





= 0, (2-1) 


where sis the direction of the flow and V ts the velocity, remains constant with inte- 
gration. With the more common Eulerian methods, a different set of fluid elements is 
evaluated at each time step, representing a different distnbution of the advective pa- 
rameters. Fjortoft’s goal was to develop a scheme in which a large time increment, Ag, 
could be used. In regions of strong wind shear or after long integration periods, however, 
the mitual grid can become greatly distorted, as illustrated in Fig. 1. Data points may 


cluster in a relatively small area, leaving large areas unanalyzed. 


> 


Figure 1. Distortion of Originally Square Grid with Time 


To avoid the problem of distortion and still make explicit use of equation (2-1), 
Win-Nielsen (1959) proposed a quasi-Lagrangian, or trajectory, scheme. The fluid’s 
properties are evaluated on a regular, stationary grid at each time step by determining 
the off-grid positions from which the fluid elements originated. 

The semi-Lagrangian method, first proposed by Sawyer (1963), is a modified version 


of the Wiin-Nielsen scheme. In this study, a one-dimensional, three time level version 


of the semi-Lagrangian method, similar to the two-dimensional scheme described by 


Robert (1981), is used. For each parameter, Q, such that 
U oe R 
= AF ay (2 — 2) 


where R represents all the remaining terms, the left hand side is approximated at the grid 


poemt, x, , by 


O(x;,0+ At) O(x,—2a,t — Ad Fs 


2At vin@ 
Palere, OF the midpoint rule, 
a, = Ate U(x; —a,.0). (2 —4) 
The a’s must be solved for iteratively. For each n, 
ay = Al < Cea, oon). (2-5) 


Subtracting equation (2-4) from (2-5) and applving the Lipschitz condition leads to 


la; —a|| Eo. 


= Are| 





| (2 — 6) 


la" —ail 
Thus, equation (2-5) will converge whenever the right hand side of Equation (2-6) 1s less 
than one. Kuo and Williams (1989) argue that no more than three iterations are needed, 
if this condition is met, since equation (2-4) is only third order accurate with respect to 
Ar and equation (2-5) becomes order Ar more accurate with each iteration. Pudvkiewicz 
and Staniforth (1984) point out that the maximum horizontal wind shear in most 
meteorological models is of the order 10~s-!, so a time step of ArS3 hours 1s necessary 
for convergence. By contrast, the time step in most Eulerian models used in 
meteorological applications is restricted by the Courant-Fnedrichs-Lewy (CFL) stability 


criterion, 


<1, (ae) 


where C is the Courant number. The time step is limited by both the maximum wind 


* - 


speed, L”_.,, and the resolution of the model. For a given U,,,,, the higher the resolution, 


Max? 


the smaller the time step: 


B. DEVELOPMENT OF HYDRAULIC JUMPS 

In this study, the one-dimensional shallow water system of equations is used to ex- 
amine the nature of the flow of an incompressible, hydrostatic, homogeneous, inviscid 
fluid over an obstacle in both a non-rotating and rotating system. Under certain con- 
ditions, it is possible hydraulic jumps may form in one or both cases. 

1. The Non-Rotating S stem 

Houghton and Kasahara (1968) studied the non-rotating shallow water system, 

illustrated in Fig. 2, and determined the conditions which would lead to the formation 
of hvdraulic jumps. When there is no Coriolis force, 1.e. f= 0, then the one-dimensional 


shallow water equations are 





Cu Cu Ci 

er : g —— — 0 2- 

Ct a Cx +8 Cx ( 8) 
and 

pp A a gp Ae = (2-9) 

Cl x 


Where uw 1s the horizontal velocity, H 1s the depth of the fluid and 


h’ = H— hyp. (2 — 10) 


aT kL 


Ch ; : 
Zy represents the height perturbation on the free surface. 

At 1<0Q, the fluid is at rest, and the free surface is flat. At 2=0, the fluid is 
given a constant velocity of u, throughout the entire domain. At t= oo, the fluid is at 


its new steady state given by 








ur ue —_ 
5 + 8(Hys + hg) = + gh = C, (2 — 11) 
and 
u,H..=u,h = Cy. (2 — 12) 





Figure 2. The Shallow Water Model 


Eliminating //, from equation (2-11) using Equation (2-12) and introducing the 


following dimensionless parameters, 





ud hp tt 
d b 


Fr=—==, R=— and UU=— (2 —13) 
aj gh Ki : 





into the resulting equation leads to 


2 2 
da fr 1)+1=0. (2 — 14) 


2 





ye etl = 


For a given Froude number, Fr, R can be plotted as a function of U. using equation 
(2-14). A plot of F versus U., shown in Fig. 3 for Fr? = 0.1, illustrates there are three 
real roots to equation (2-14) if R< Regjrjcy,. Only one of these roots is physically 
meaningful. If R> Repiricaz, Only one solution to equation (2-14) exists, but since this 
solution has no physical meaning, a hydraulic jump is anticipated. Reg;ric4, Can be ex- 


pressed as a function of the Froude number: 


2 
Rerirican = FP — > Frd +1. (2 — 15) 


A plot of Repirica, Versus Fr, illustrated in Fig. 4 indicates the three classes of solution 
to eqyuauon (2-14) idenufied by Tloughton and Kasahara. 

Domain | 1s the’subenitical range iy which K = Keo, ANG? = 4 | eectedcn 
state free surface height of a fluid flow Which meets these criteria will dip over the ob- 
stacle. The velocity will increase over the obstacle but will remain less than u., the speed 
that corresponds to the condition R = Repjrjc4,- Domain If is the ‘jump region’ in which 
R> Rerica, and a hydrauhe jump forms. Domain HI is the supercritical range in which 
R= Rearjcyz, DUt Fr> 1. At its steady stute, the iree surface height rises over tilemers 


stacle, anc the velocity decreases but remains greater than w.. 


ee CRITICAL 


a0) 


Figure 3ameK <p een Versus. U. lon) 70m! 
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2. The Effect of Rotation 


If the coordinate system is rotating, 1.e. f#0, then the shallow water system 


equations for a flat lower boundary and no y dependence are 


Gil Ou OH > » 
io ace a (2-08) 
RON coe Zs > 

Eps a + fu=0 (2 — 17) 


and 


Eu =o (2 — 18) 
Cl Cx CX 


If the variables are scaled as follows: 


He 





























x= 1s, (= Ue, 
V Bh, 
u=Uu, ve(42)y, He R+Us/4 a, (2 — 19) 
oO 
a etal ae — == 
i gh 
then equations (2-16), (2-17) and (2-18) become 
‘ P ; 
ale 3 ee (2 =a 
Ol; ONS R 
. 
et ae (2 — 21) 
and 
Cie Ch, Cu, Cite > oe 
ae + Flu, EFS +h, oe )+ om = 0), (2 — 22) 


Unless the Coriolis term dominates, hydraulic jumps may form due to the nonlinear 


advection terms. The no-jump condition can be expressed as 





Or 
F=A+R>, (2 =24) 


where 4 is some number much larger than one. The Williams and Hori (1970) investi- 


gation found that A is in the range of 6.5 to 7.0. 


3. Scale Collapse Problem 
Computational difficulties may arise When a numerical model must represent a 
physical process in which a large gradient or discontinuity forms. Kuo and Williams 
(1989) considered a simple scale-collapse case, for which the governing equation was the 


one dimensional nonlinear advection equation 


ft 4 Sh 9, (=) 
Cl Co 
and the initial condition 
u(x), 9 = — tan” (x — x): (2 — 26) 


This problem can be solved for analytically. and that solution leads to 


Cu l - 

——.> cine oc oas fl... ee) 

Cx a, a 7s | 

Kuo and Willams compared the second-order centered finite difference, the 
Chebvshev-Tau and the semi-Lagrangian solutions of equation (2-25) and (2-26). The 
scale-collapse at x = x, produces large, fluctuating errors throughout the entire domain 
when either the finite difference or Chebyshev-Tau method is used. The sem- 


Lagrangian method responds differently. As ¢ approaches one from below, 
esi (2 — 28) 
Gx 


for iAx=x,. Therefore, equation (2-5) will not converge in this region. However, if a 
fixed number of iterations (m = 3) 1s used to solve for the a's, the errors which result 
when and where equation (2-28) is true remain near the scale collapse zone. Thus, the 
semi-Lagrangian method 1s still useful even after discontinuities develop. 

In this study, two semi-Lagrangian schemes for the shallow water equations are 
compared with the traditional leap-frog finite difference method. In the first, the R terms 
in equation (2-2) are treated explicitly; that is, they are evaluated at (x,—a,v). This 
scheme is expected to handle large shear zones, but the freedom from the CFL criterion 
cannot be exploited, because the time step is still limited by the gravity wave terms. The 
second scheme 1s a semi-Lagrangian, semi-implicit scheme similar to the one presented 
by Pudykiewicz and Staniforth (1984), but modified to include topographical effects. 


The terms responsible for the gravity waves are treated implicitly by averaging between 


(x,¢+ Az) and (x,—2a,.—An). The remaining terms are evaluated at (x,—a,,t). The 
semi-Lagrangian, semi-implicit scheme requires more computational effort, but it per- 


mits a much greater time step. 
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IY. FORMULATION OF THE MODELS 
A. THE GOVERNING EQUATIONS 


The governing equations for the shallow water model, illustrated in Fig. 2, are 











Cll Cu CH Che 
Btw a of ax Ax )-fr=0, (3 — 1) 
ey ey CH Che 
Bu a Ris By )+fem0 (37 2) 
and 
ral. F Oi re “aie 
A ~ ax (ull) + ear eh) =O, (33) 


In the one-dimensional model, the mean flow, F’, 1s zonal and in geostrophic balance 
with the mean free surface height, such that 


Oo 
u — fees rR 3-4 
: cy ( 
and 
g ch ° 
y= — — = Q), 3-5 


and perturbations in the dependent variables, u’, v’ and A’ are in the x direction only. 


Multiplving equation (3-3) by g and substituting 


u(x,t) =u, + u'(x,0), C=s) 
u(x,t) = v'(x,0) (3-7) 

and 
$(x,1) = 60) — ba(xy) + "(2,0 = 8H (3 — 8) 


into equations (3-1), (3-2) and (3-3) vields the following set of equations 


It 

















se +(u, +0) So + = -fv' = Oz 
aaa +1) ae (One!) 
and 
Stu, + uh( = a ae )+@-6,+ 0 - )+ 
{ B-Se } -0 @e ui 
Cy cy 


For convenience, the bottom topography 1s assumed to be parallel to the mean free 


surface height in the y direction, so that 














—— _ (3 — 12) 
oT Cy 
and equation (3-11) reduces to 
Cp’ C Ons EGG ps) ye Cu 
5 +(u,+u)l aaa ee re ao) — = 0 3 —ids 
: (td, ) eS ane One ) Are ( ) 


B. FINITE DIFFERENCE SCHEME 
One of the simplest numerical methods is finite differencing centered in time and 
space, the leapfrog scheme. On the staggered grid, equations (3-9), (3-10) and (3-13) 


approximated by this method become 


u'(x.t+ At)— "(x0 — Ad) u(x + Axa) — u'(x — Ax. 
7, Se = x0)) Ax = 
Meee ye Hie Big ple ES 
eee ee Hf es oo = G14 
Ax p 


w(x + Ar) — v'(x,t—- Ad ie also (208 (0) eos sf) 
2At aes) 2AN = 


u(x + Ax N+u(x- ass f) 
- a (Srl) 





~ 


and 


O'(x~+tAnd—P'(xt—AdD _ p'(x + Axa) — b'(x — Ax.n 
2a a) PAX 


7 O p(x + Ax) — bp(x — Ax) ; 
DI 


2Ax 
_ u'(xt ge Yo—u(x- oe f) 
— (6 - dg(x) + i —— | Bec) 
where 
u(x — 2S tulle t Ax BE) 
fi(x.t) =u, + _ | — i 
we 


C. SEMI-LAGRANGIAN, EXPLICIT SCHEME 

The second numerical scheme uses the semi-Lagrangian method described bv Robert 
(1981), modified for one dimensional flow over a surface with topography. The re- 
maining terms are centered in space about (x — a,z). Equations (3-9), (3-10) and (3-13) 
approximated by this method become 


u(x. + At) — u(x — 2a— Ad) Og’ _ 
a Ro, Ct 


v(x + At) — v(x — 2a,t— Ad a 


sar — ful(x — a,t) (3 — 19) 


and 


b'(x.t+ At) — $'(x — 2a,t— Ad) @ p(x) — Op(x — 2a) a 
2Ar 7 2Ar 7 
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a (6 aa Pp = op’) Se leven) e = 20) 








where 
A. Be 
SOCOM ot a ee 
wo = (3 — 21) 
Equation (3-21) is used to estimate = and ae at the grid points, and then an in- 


terpolation scheme is used to evaluate these terms at (x—a). An alternative m -hod 
would be to approximate the derivatives immediately at (x —a) with the following 


equation: 


Although this method appears more direct than using equation (3-21), it requires twice 
as manv interpolations and therefore more computational effort. 


Note, 


COgp COZ Cbg — Og(x) — Ga(x — 2a) 


Gx Ct Cx QAI 











| 

at 
a | 
2 


(3 — 23) 


“~ 


, | . . Sere 
since there 1s no local time change in the bottom topography, 1.e., ay 1S Zero. 





D. SENII-LAGRANGIAN, SEMI-IMPLICIT SCHEME 

This scheme is similar to the semi-Lagrangian, explicit scheme presented above ex- 
cept the gravity wave terms are treated semi-implicitly. The development of this scheme 
is similar to the development of the semi-Lagrangian, semi-implicit scheme developed 
by Staniforth and Temperton (1986) for a flat bottom system. The finite difference 
equations now become 


u'(x,0+ At) — u'(x — 2a,t — Ad) 1 | é¢’ 6d’ 
aa a — =a — 0 loxtA0 ae a aaa 
+ fv'(x —a,0), (3 — 24) 


v'(x,0+ Ar) — w(x — 2a, — Ad) 


ie = —fu'(x — a,t) (3 — 25) 
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and 


O'(xi+ Ad — o'(x — 2art— Ar P plx) — Op(x — 2a) d 7 
2Atr 7 2At ee AG I(x—a,2) 





=F [GF - b5(2)) BE harean + 8 ~ dsl - 24) HE I reseay | B28) 


The divergence term at (t+ Az) 1s eliminated from equation (3-26) by first finite differ- 


entiating equation (3-24) with respect to x, which vields, 


2 z ’ rae 
ae aus ae ale ee 
Ox (x.t+A0) Ox (x—2a,t—Ar) 5x2 (x,t+An 5x2 (x—2a,f—At) 


_ 


+ 2Aif-— easy (3 = 27) 





where 


56’ bx t+ Ax) — 26'(x) + o'r — Ax) 


ae? 2 : 
Ox Ax 





(32) 


then substituting equation (3-27) into (3-26) and collecting all (x,:+ Az) terms on the 
left-hand side 


a 5’ Le on A =_> — oe a. ee ’ _ > f= A 
( ix? APG — bx) ) nie. ( bx? APG — bax) ) iil 


_ bale) = a(x — 2a) ( A, $= Oa(x = 24) 3 bul 





bt (eee 1—Ar) 


APG — b(x)) Ar b — b,(2) 
26'(x — a) bu’ pe 
Aud — bs(x)) OX ena = Phisnas ee 


Or, expanding, 


( Ax? 
ao) 
7 Ar(S — d(x) 


Ax? 


Joven + @'(x — Ax) 


eky: = 


Ax? | 
‘(x — 2a + Ax) =| 2 = SS Soe IO OS 2 eee 
p'(x ise A) ( Ap oan om Jere a)+ $'(x — 2a x) 


D-II = iio! "Cll “LEliuaaca. (cc. 
Ax? 





a ; A. ‘ : 

b a(x) — box — 2a) ( 1 $— p(x — 2a) | u(x ~ 2a +S) — u(x 2a) 

Tw | <<< 
Ar(d — bp(x)) Ar And — (x) Ax 





2h'(x — a) u(x ~a +S) —ul(e—a- 
| Se , 
Ar(p — (x) Ox 


vixtat Ax)—v(x-—a— Ax) P 
+7) eee | F (3 — 30) 
Equation (3-30) can be written in matrix form as 
LA}{O';} = (B}, Cir). 


where {B,} are all the term: on the right-hand side of equation (3-28) and the matrix A 


iS 


O19 ell Oc 0 
lL ‘Goel. a0 0 
0 or 0 ee SOs Se 
0 Ue 
where 
2 
un (245th | (3 — 33) 

Ar'(o — @p(iAx)) 


To solve equation (3-31), it is necessary to invert the matrix A : 


(= ie te (3 — 34) 


However, since matrix A has no time dependent elements, the inverse, [A]-', must be 


solved for only once. 


E. THE SPATIAL GRID 

u’, v’ and @’, are evaluated on the staggered grid, scheme C of Arakawa and Lamb 
(1977), illustrated in Fig. 5. The use of scheme C further limits the maximum time step 
for the explicit models, already restricted by the CFL criterion, by a factor of two, as 
Haltiner and Williams (1980) demonstrated with the leapfrog finite difference scheme. 
However, they also showed that scheme C simulates ‘he shorter waves far better than 
the unstaggered grid, an important feature when investigating orographic effects. Cyclic 
boundary conditions are used to avoid the difficulties associated with the advection of 
fluid elements from outside the domain onto the grid. 

‘An interpolation scheme is needed to evaluate the dependent variables at (x — a) and 
(x — 2a). but the resultant smoothing of the variables between grid points can represent 
significant errors. Kuo and Wilhams (1989) have shown that in regions where the 
Courant number 1s less than one. interpolation errors dominate time truncation errors. 
In these situations, the forecasts become miore accurate as Ar is increased, because fewer 
time steps require fewer interpolations. 

In this study, a cubic spline interpolator is used. Although computauonallv expen- 
sive, the cubic spline, because of it accuracy, has been a popular choice for use with the 
semi-Lagrangian method (Robert (1981), Pudykiewicz and Staniforth (1984), Ritchie 
(1987), Kuo and Williams (1989).) The global nature of the cubic spline ( all grid values 
are used to estimate to a single point), led Pudykiewicz and Staniforth (1984) to suggest 
that this scheme would be inappropriate where locally steep gradients occur and to 
speculate that a local interpolating function should be used in these situations to mini- 
mize the Gibbs’ phenomenon. However, Kuo’s and Williams’ (1989) investigation of a 
scale collapse model found that the deviations in semi-Lagrangian forecast from the 
analytical solution remained small and confined to the region of the large wind conver- 


gence, even though a cubic spline interpolator was used. 


NT 


(2,v) | 





igure 5. Staggered Grid 


lV oon i~P PRIMIENTSe¥ND RESULTS 


In the first experiment, the three models’ forecasts for flow over a topographical 
ridge in a non-rotating system are compared. The same cases that Petroliagis (1988) 
investigated with a Galerkin model are examined here. All three domains identified by 
lfoughton and Kasahara (1968) are represented. 

In the second experiment, the effect of rotation upon topographically induced jump 
formations ts examined. The three iodels are run with the same sets of parameters used 
in the first experiment, except /= 0.000103 sec-', corresponding to the Coriolis force at 
45° N. 


A. INITIAL CONDITIONS 
The domain of integration, L, 1s 7276 km in the x direction, divided into 48 incre- 
ments of length aS wis evaluated on the odd numbered grid points; v and @ are eval- 


uated on the even numbered grid points. The boundary conditions are cyclic, sueh that 
2(O) = u(L), (0) ae) andua (0) — (1). (5 — 1) 


The ridge profile, illustrated in Fig. 6 is described by 


Ri 20 
Hp sin“ = for S20 Ae = 324% 


0 elsewhere. 





Figure 6. Profile of the North-South ridge 
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The five cases examined, summarized in Table 1, are parameterized by the mean 
height, # , the ratio of the mountain peak to the mean height, R, and the Froude num- 
ber, Fr. These are the same sets of parameters that Petrolagis (1988) used. However, 
for these experiments L and Ax are increased tenfold which increases the time scale, so 


that f could be important in the second set of experiment. 


Table 1. PARAMETERS FOR CASES I THROUGH V 


iim [| Rk | | uma] D1 
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A five-minute time step, A/;y, is used with the Finite Difference Explicit (FDEX) and 
Semi-Lagrangian Explicit (SLEX) models. A one-hour time step. Av,,, 1s used with the 
Semi-Lagrangian Semi-Implicit (SLSI]) model. The Semi-Lagrangian Semi-Implicit 
model is also run at a higher resolution (HRSL) for comparison with the other models. 
For HRSL, Ax is reduced by a factor of four, and Avs,, is used. 

Atry 1s limited by the CFL stability criterion based on the speed of the explicitly 


treated external gravity waves. The mean speed of these waves, c, is given by 


=|(Frt vo |. (4-3) 
A mean Courant number for the explicit cases, Czy can be defined as 


cAtry 


2 TAs 








As discussed in chapter 3, the factor of two arises because scheme C is used. Cry is ap- 


proximately 0.25-0.30 for the cases considered. 


B. RESULTS I 
The wu field for case I after two, four, six and eight hours of integration are shown 


in Fig. 7. All the models produce a stationary speed maximum centered over the ridge, 
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Which the theory predicts for the subcritical range (Domain I). FDEX and SLEX 
slightly underdevelop the speed maximum while SLSI is in good agreement with HRSL. 
All the models show the development of a secondary maximum that propagates east- 
ward. After eight hours the downstream disturbances have propagated around to the 
east side of the ridge and begin to interfere with the upstream fields. 

As anticipated, a low pressure region develops on the lee side of the mountain 
(Fig. 8). The high pressure field to the east of the ridge remains fairly stationary in the 
SLSI and HRSL forecasts. This feature drifts upstream perceptibly in the FDEX and 
SLEX forecasts Theory predicts the formation of a wave train on the lee side of the 
mountain. FDEX, SLEX and SLSI all lag behind the high resolution model in devel- 
oping these waves. SLSI is only slightly behind. The FDEX becomes appreciably 
damped, but SLEX has the poorest results due to excessive smoothing after numerous 
interpolations. 

Case II 1s also in Domain I, however the mean flow, Uy;z,,, is stronger and the 
mountain peak, //;. 1s higher. Although the perturbations are relatively larger. the u and 
p’ fields, shown in Fig. 9 and Fig. 10 for two, four, six and eight hours, are quite similar 
to their counterparts in case I. 

Case III and case IV le in Domain II, so hydraulic jumps are anticipated. Fig. 11 
and Fig. 12 illustrate the u and @’ fields at three, six, nine and twelve hours for case III. 
All models forecast a wind speed maximum over the ndge. The FDEX solution does 
not appear in Fig. 11C and Fig. 11D and in Fig. 12C and Fig. 12D, because FDEX is 
unstable, failing to converge after six hours, due to the large nonlinear interactions in the 
vicinity where the jump is developing. The HRSL model continues to deepen and 
broaden the pressure minimum on the lee side of the ridge until, after eleven hours of 
integration, the pressure has dropped to zero. Once the free surface height hits the bot- 
tom topography, the integration is stopped. The perturbation fields develop slightly 
slower in the SLSI model; the free surface hits the bottom after 13 hours of integration. 
Initially, the SLEX forecasts are very similar to SLSI’s. With time, however, the SLEX 
perturbations dampen out, again due to interpolation error. 

The u and @’ fields at three, six, nine and twelve hours for case IV are shown in 
Fig. 13 and Fig. 14. The amplitude of disturbances are increasing, but because of the 
cyclic boundary conditions, there is interference from the transient part of the solution 
before a jump can form. 

Case V is in Domain III, the supercritical range. The u and ¢’ fields from all the 


models are compared in Fig. 15 and Fig. 16. SLSI again agrees well with HRSL. Both 
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produce the speed minimum over the ridge as the theory predicts. FDEX and SLEX 
develop this minimum slightly upstream. All the models show a pressure minimum to 
the lee of the obstacle. The SLSI and HRSL remain rather stable, until the transient 
solutions propagating to the east circle back around to the west of the ridge. The SLEX 
model exhibits exceesive smoothing from accumulating interpolation errors, as it does 
for all the other cases. The FDEX solution becomes unstable in this case, unable to 


handle the nonlinear interactions when Fr exceeds one. 


C. RESULTS Il 

W en the system is allowed to rotate, perturbations arise in v as well as in u and @. 
These disturbances in v are out of phase with the uw’ field, as illustrated in Fig. 17 for 
case I. After several hours of integration, perturbations in u are noticably smaller and 
the pressure drop on the lee side of the obstacle appears not as deep when compared to 
the non-rotating cases, as in Fig. I8A and Fig. 18B for the HRSL results for case If. 
Although the Coriolis force is deflecting some energy into the y direction. cases II] and 
IV are sull unstable. For both cases, the amplitudes of the disturbances in the vicinity 
of the ridge continually increase. This is illustrated in Fig. 19 for the total wu field in case 
IV after three, six, nine and twelve hours of integration. 

Fig. 19 1s also representative of the relative performance of the models for all cases 
in the rotating system. Not suprisingly, the results of these comparisons are much the 
Same as in the non-rotating cases. SLSI agrees most closely with the high resolution 
model. FDEX 1s unable to capture all the detuils in the disturbance fields and has a 
tendency to place the extrema slightly upstream from the HRSL model. The SLEX 


forecasts are the poorest due to the dominance of interpolation errors. 
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Figure 7. Experiment 1, Case I: Total u Field (sn/s) after A) Two B) Four C) Six 
D) Eight Hours of Integration for HRSL (solid line), SLSI (dashed line), SLEX 
(broken line) and FDEX (dotted line). 
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Figure 8. Experiment 1, Case I: @’ Field (s/s)? after A) Two B) Four C) Six D) 
Eight Hours of Integration for IIRSL (solid line), SLSI (dashed line), SLEX (broken 
line) and FDEX (dotted line). 
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Figure 9. Experiment I, Case IJ: Total u Field (s/s) after A) Two B) Four C) Six 
D) Eight Hours of Integration for HRSL (solid line), SLSI (dashed line), SLEX 
(broken line) and FDEX (dotted line). 
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Figure 10. Experiment I, Case II: @’ Field (m/s)? after A) Two B) Four C) Six 
D) Eight Hours of Integration for HRSL (solid line), SLSI (dashed line), SLEX 
(broken line) and FDEX (dotted line). 
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Figure 11. Experiment 1, Case III: Total u Field (m/s) after A) Two B) Four C) 
Six D) Eight Hours of Integration for FIRSL (solid line), SLSI (dashed line), SLEX 
(broken line) and FDEX (dotted line). 
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Figure 12. Experiment 1, Case II]: @’ Field (s/s)? after A) Two B) Four C) Six 
D) Eight Hours of Integration for HRSL (solid line), SLSI (dashed line), SLEX 
(broken line) and FDEX (dotted line). 
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Figure 13. Experiment 1, Case 1V: Total u Field (m/s) after A) Two B) Four C) 
Six D) Eight Hours of Integration for HIRSL (solid line), SLSI (dashed line), SLEX 
(broken line) and FDEX (dotted line). 
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iment I, Case IV: @’ Field (m/s)? after A) Two B) Four C) Six 
f Integration for HRSL (solid line), SLS] 


(dashed line), SLEX 
/ DEX (dotted line). 
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Figure 15. Experiment 1, Case V: Total u Field (s/s) after A) Two B) Four C) 
Six D) Eight Hours of Integration for EIRSL (solid line), SLSI (dashed line), SLEX 
(broken line) and FDEX (dotted line). 
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Figure 16. Experiment 1, Case V: $’ Field (m/s)? after A) Two B) Four C) Six 
D) Eight Hours of Integration for HRSL (solid line), SLSI (dashed line), SLEX 
(broken line) and FDEX (dotted line). 
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Figure 17. Experiment 2, Case I: uw’ and v’ (s/s) after A) Four B) Eight Hours of 


Integration with TIRSL. 
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Figure 18. Comparison of Perturbations in Rotating and Non-rotating System (#1/s)? 
for Case IT A) wand B) f’ with HRSL. Solid line indicates non-rotating system; the 


dashed line indicates soluuion for the rotating system. 
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Figure 19. Experiment 2, Case IV: Total u Field (sn/s) after A) Two B) Four C) 
Six [ID) Eight Hours of Integration for HRSL (solid line), SLSI (dashed line), SLEX 
(broken line) and FDEX (dotted line). 
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V. CONCLUSIONS 


This investigation is a direct comparison of a senu-Lagrangian, semi-implicit model 
to the semi-Lagrangian, explicit and leapfrog finite difference models for the shallow 
water system with bottom topography. Of particular interest is the models’ ability to 
simulate hvdrolaulic jumps. 

In the first experiment, the Coriolis parameter is set to zero, and the cases which 
Petroliagis (1988) exanuned with a Galerkin finite element model are reconstructed, al- 
beit on a larger domain. Initially, the flow was uniform and zonal, and the free surface 
was flat. -\ecording to the theory of Houghton and Kasahara (1968), the first two cases 
are subcritical. A stable speed maximum should appear over the ridge. The last case 1s 
supercritical, and a speed minimum should develop over the ridge. The remaining cases 
are unstable and should produce a jump in the vicinity of the obstacle. 

The SLSI forecasts are superior to the [DEX and SLIEN foreeasts iorealineas. 
SLEN forecasts are destroved by excessive smoothing due to interpolations; FDEX 
breaks down when the non-linear interaction become too large. SLSI results. however, 
are consistent with hydraulic jump theory. In the no junnp cases, it forecasts stable speed 
and pressure maxima and minima where thev are predicted bv theory. SLS1 forecasts 
increasing Winds and dropping pressure for the jump cases. !flowever, the details of the 
hvdraulic jump which forms in case III are only apparent when the resolution 1s 1n- 
creased, and in Case 1V, the domain is too small for a hydraulic jump to form before 
interference from the transient solutions due to the evelic boundary conditions occurs. 
Further testing should be done on a larger domain (more tsolated ridge) with the higher 
resolution semi-Lagrangian, semi-implicit model, IIIRSL. 

In the second experiment, the same cases were re-examined, but the Coriolis pa- 
rameter for 45°N was used to examine the effect of rotation on jump formation. Clearly, 
rotation slightly suppresses the perturbations in the free surface and delays the f sma- 


tion of hydraulic jumps, but further study is necessary for a more quantitative analysis. 
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